Surface defects and conduction in polar oxide heterostructures 
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The polar interface between LaA103 and SrTiOs has shown promise as a field effect transistor, 
with reduced (nanoscale) feature sizes and potentially added functionality over conventional semi- 
conductor systems. However, the mobility of the interfacial two-dimensional electron gas (2DEG) 
is lower than desirable. Therefore to progress, the highly debated origin of the 2DEG must be 
understood. Here we present a case for surface redox reactions as the origin of the 2DEG, in partic- 
ular surface O vacancies, using a model supported by first principles calculations that describes the 
redox formation. In agreement with recent spectroscopic and transport measurements, we predict 
a stabilization of such redox processes (and hence Ti 3d occupation) with film thickness beyond a 
critical value, which can be smaller than the critical thickness for 2D electronic conduction, since 
the surface defects generate trapping potentials that will affect the interface electron mobility. Sev- 
eral other recent experimental results, such as lack of core level broadening and shifts, find natural 
explanation. Pristine systems will likely require changed growth conditions or modified materials 
with a higher vacancy free energy. 



I. INTRODUCTION 

Complex oxides offer the potential to replace con- 
ventional semiconductors in a range of devices due to 
re duce d feature sizes and added functionality (see e.g. 
ref . The polar interface between LaA10 3 (LAO) and 
SrTi0 3 (STO)^has shown promise as a field effect tran- 
sistor SIS. One problem hindering its development is the 
low mobility of the interface two-dimensional electron 
gas (2DEG). To progress, the origin of the 2DEG must 
be understood. Explanations proposed to date can be 
classed into three categories: (i) electron transfer coun- 
tering the "polar catastrophe'^, (ii) doping through O 
vacancies in the filnPEHU or j n the su bstrate^ H12l, and 
(in) cation intermixing at the interface 1 6 * 13 * 14 !. Growth 
conditions have been shown to affect the observed behav- 
ior, particularly O2 partial pressure and post-annealing 
treatments, affecting carrier density and confinement to 
the interf ace ( as opposed to conduction through the STO 
substrate) 4 ^. 

The polar catastrophe^ arises from a polarization dis- 
continuity between the non-polar STO substrate and po- 
lar LAO film's] i n LAO films grown on Ti02 terminated 
STO substrates this polar discontinuity induces a diver- 
gence of the electric displacement field (V • D) in the 
pristine system equal to considering the effective polar- 
ization response and net effective charges a c of precisely 
+e/2 and — e/2 per formula unit at the interface and sur- 
face, respectively^* 1 ^ 1 . The electrostatic potential then 
builds up across the LAO film, accumulating electrostatic 
energy. To counter it a charge transfer between the inter- 
face and the surface is required. A mechanism intrinsic 
to the pristine system is given by the transfer of electrons 
from the surface to the interface. From band bending ar- 
guments, electrons are transferred from the O 2p at the 
top of the valence band at the LAO surface, to the Ti 



3d conduction band at the interface, once the potential 
drop across the LAO layer reaches the effective band gap, 
which is calculated to happen for a LAO thickness of five 
unit cells (see ref and within) . However the absence 
of a 2D hole gas at the surface and the observation of 
populated Ti 3d states for films as thin as one or two 
bilayers^Hll] raises doubts about this mechanism. 

Two important points should be kept in mind. Firstly, 
the samples are generally not kept or measured in vac- 
uum. The surface can thus be covered by adsorbants like 
water, and it will likely be far from an ideal surface ter- 
mination, possibly including chemical alterations such as 
the hydroxylation seen on many wet oxides. The second 
point is more central, however: whatever the chemistry, 
the relevant electrostatics across the film can only be af- 
fected by processes in which charge is altered at each side 
of the film, either by charge transport across it, or by 
charges arriving at either side from external reservoirs. 
If there are no such external sources and the chemical 
processes are confined to the surface, the remaining pos- 
sibility is that of surface redox processes. They trans- 
form surface bound charge into free-carrier charge, the 
electrons or holes then being free to move to the buried 
interface. The clearest and quite relevant example is that 
of oxygen vacancy formation whereby surface 2 ~ an- 
ions transform into O2 molecules, releasin g tw o electrons 
to the n in terfac e, as illustrated in Fig. lEEHU Surface 
protonatior l 22 * 23 * is an analogous process in which the sur- 
face is also reduced by the oxidation of O 2- into O2, al- 
though in this case the process depends on the presence of 
water, H 2 0-s- 1/2 2 + 2e~ +2H+, and both the energet- 
ics and kinetics will be different from the previous one^l. 
The protons attach to the surface O atoms, while the 
electrons are again free to go to the n interface (note that 
the non- redox hydroxylation, H2O— > OH~+H + , does not 
affect the electrostatics across the film). 
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In section II we present a model for the formation of 
surface redox processes. We then apply the model to the 
LAO/STO system in section III using parameters ob- 
tained from first-principles calculations (see appendix). 
We show that, for such a surface redox process, (i) the 
process is favored by film thickness through the effect 
of the electric field across the film; (ii) the growth of 
the density of the related surface defects with film thick- 
ness is predicted, showing a minimum critical thickness; 
(Hi) the model is in excellent agreement with first prin- 
ciples calculations and reasonable agreement with what 
is observed experimentally for Ti 3d occupation; (iv) the 
potential drop across the film does not change substan- 
tially with thickness once vacancies start to appear; (v) 
for thin films the carriers at the interface are still trapped 
by the electrostatic potential generated by the vacancies 
in deep double-donor states: a strongly disordered two- 
dimensional electron system; (vi) as thickness increases 
the levels become shallower and the number density in- 
creases, closing the gap to the conduction band; (vii) the 
conduction onset is therefore at a higher film thickness 
than the one for interfacial carrier population. Finally 
we note that the model can also be applied to p-type 
interfaces and ultra thin ferroelectric films, where alter- 
ation of the surface chemistry has also been proposed as 
a possible screening mechanism (see for example ref .^) . 



II. MODEL 

We consider a pristine polar thin film on a non-polar 
substrate defining a n-type interface (see Fig. 1), and 
a surface reduction process (the model can be trivially 
extended to p-type interfaces and surface oxidizing re- 
actions, see below, and to ferroelectric thin films). The 
introduction of a surface defect via a redox reaction pro- 
duces a donor level in the gap. The defect provides Z 
electrons (two for an O vacancy) that can transfer to the 
interface (see Fig. 1), thus contributing to the screening 
of the polarization or compositional charge a c at either 
side of the polar film (±0.5e/f.u. in LAO^^. 

Considering this electron migration, we can model the 
formation energy of one such surface defect, Ef, in the 
presence of an area density n of surface defects as 



(a) 



E f (n) = C + E e (ri) + an, 



(1) 



where we have separated an electrostatic term associated 
with the internal electric field in the film, Eg(n), from a 
surface/interface chemistry term, C, and a term account- 
ing for defect-defect interactions other than electrostatic 
in a mean-field sense. It can be seen as arising from 

others 

E f (n) = E f -Z(W-E CD ) + E e (n)+ An, ( 2 ) 



where E® is the formation energy of an isolated surface 
defect in the absence of a field across the film, and where 
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FIG. 1: (Color online) Schematic band diagram of an inter- 
face between a polar film and a non-polar substrate along the 
normal direction, z. (a) The pristine system under the crit- 
ical film thickness, (b) The creation of a donor state at the 
surface via a redox reaction and subsequent electron transfer, 
(c) The reconstruction reduces the film's electric field. 



W and E CD are defined in Fig. ]P. E° f (and thus C) 
depends on the particular surface chemical process, and 
the reference chemical potential for the relevant redox 
counterpart species in the environment, e.g. fxo 2 (P,T), 
which depends on experimental conditions. 

Taking Eq. 1, the surface excess energy for a given 
area density of surface defects is then 

r n i 

Q(n) = / E f (n')dn' = Cn + fl e (n) + -an 2 . (3) 
Jo 2 

The key of the proposed energy decomposition is that the 
Of term is simply the energy gain of partly discharging 
a capacitor, which is 



d 

2e 



(4) 



where d is the film thickness, e is the LAO dielectric con- 
stant, a c is the compositional charge, and a v = n{Ze) is 
the charge density of the carriers confined to the interface 
(note that these electrons may not all be mobile, as dis- 
cussed below). Eq. 4 assumes no screening by electronic 
reconstruction, which is right if the onset for defect stabil- 
isation happens earlier than the one for electronic recon- 
struction. A complete description of all possible regimes 
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will be presented elsewhere^. We limit ourselves to the 
regime given by Eq. 4 since a wider discussion of the 
model is irrelevant here. 

The equilibrium defect density is determined by finding 
the value that minimizes Q. Taking Eqs. 3 and 4, 

dZe<j c — Ce 
(Ze) z d + ae 

A critical thickness arises for defect stabilization, 

<k = Ce/(Zea c ), (6) 

n tending to a c /Ze for large d, which is the value required 
to completely screen the film's intrinsic polarization. 

III. DISCUSSION 

A. Interface carrier density 

We now consider the specific case of O vacancy for- 
mation at th e sur face as the most prominent candidate 
redox process^^l. Fig 2a) shows quantitative agreement 
between the model's Ef = -O(n) = i L Efdn' and 
first principles calculations of the surface vacancy forma- 
tion energies in refP^on the full LAO/STO structure. Et 
is the right magnitude to compare to first-principles re- 
sults since it accounts for the energy difference between 
the system with a given concentration of surface defects 
(n) and the pristine system (n — 0), per surface defect, 
thus Q(n)/n. The physical constants used in the model 
were determined independently by separate DFT calcu- 
lations (see appendix) and then compared to DFT results 
for films of varying thickness (Fig. 2a)). The predicted 
behavior of n (and a v ) in LAO/STO is shown in Fig. 
2b), where it is compared with Ti 3d occupation (both 
trapped and mobile) as measured with HAXPESP^. Two 
bands are plotted, the colored one uses a range in C, 
the striped one a range in a and e (see appendix) . The 
bulk dielectric constant of LAO is about 24^. It may be 
substantially different for a strained ultra-thin film (see 
refpSl and within) and so the range 21 < e < 46 has been 
considered. 

The agreement in Fig. 2b) between model and experi- 
ment is only qualitative given the ambiguities in some of 
the magnitudes of key parameters defining the problem, 
most notably the chemical potential of O2 in experimen- 
tal conditions. Despite this, the model predicts a critical 
thickness for the appearance of carriers at the interface 
for a LAO film thickness below the five unit cells pre- 
dicted by the purely electronic mechanism. Other quali- 
tative features observed but not understood in this sys- 
tem also find a natural explanation (below). 



B. Electric field in LAO: pinning of potential drop 

The electrostatic potential drop across the LAO film is 
V=(a c - a v )d/e. (7) 
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FIG. 2: (Color online) a) Defect formation energy, Ef (see 
definition in text) versus LAO film thickness d for various va- 
cancy densities n. The model (lines) is compared with the 
DFT calculations (circles) of refP of the surface vacancy for- 
mation energy on the full LAO/STO structure (see appendix 
for the determination of the model parameters) b) Equilib- 
rium area density of interface carriers a v versus d. The red 
(grey) band is the model prediction for 2.1 eV < C < 5.0 
eV, e = 25 and a = 0.8 eV/(vac/f.u.) 2 . The striped band is 
for C = 3.6 eV, 21 < e < 46 and < a < 8 eV/(vac/f.u.) 2 . 
The circles indicate the Ti 3d occupation as measured with 
HAXPES in refP^. Open circle indicates the sample was not 
annealed. The crosses indicate the carrier density from Hall 
measurements in refP. 



Substituting a v , the drop is essentially independent of 
thickness, V ~ C/(Ze), when the vacancy- vacancy in- 
teraction is small, a <C (Ze) 2 d/e. Using the parame- 
ters for LAO/STO the difference in potential drop per 
LAO layer added is between 0.0 and 0.2 eV/f.u., much 
smaller than the predicted thickness dependence for the 
electronic screening model, and consistent with recent 
XPS measurements, which show no core level broaden- 
ing with film thickness^. The pinning of V is also con- 
sistent with reduced cation-anion relative displacement 
with increasing LAO thickness as measured by SXRTJsH. 
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C. Onset of conduction: electron trapping 

The redox processes proposed above explain the ab- 
sence of hole-mediated transport at the surface, while 
electrons allow 2D conduction at the interface. An im- 
portant observation that remains unexplained, however, 
is the fact that the onset of interfacial Ti 3d occupation, 
as measured with HAXPESJ^HS] happens at lower film 
thickness than the onset for interfacial 2D conduction^. 
It has been suggested that the 2DEG lies in several Ti 3c? 
sub-bands, some of which are not mobile due to Anderson 
localization 29 . Whether Anderson localization occurs on 
an energy scale as high as room temperature depends on 
the energy scale of the disorder distribution. The surface 
defects associated to the redox processes represent point 
sources of effective charge, very much as a dopant in a 
semiconductor^, e.g. +2 e for an O vacancy. They then 
generate trapping potentials for the carriers at the inter- 
face plane of the form Vtrap = Ze 2 /e\/p 2 + d 2 , in atomic 
units, where d is the film thickness and p 2 = x 2 + y 2 
corresponds to the radial variable in the plane. This po- 
tential is sketched in Fig. 3 for several d values. Its depth 
decays with thickness as 1/d. Fig. 3 shows estimates of 
the ground state electron level associated to the double 
donor state arising at the interface due to an O vacancy 
at the surface. These trapped interface levels may be the 
'in-gap states' seen in a recent spectroscopic study^. For 
a thin film the traps are deep and few, but as it grows 
thicker, the donor states become shallower and the area 
density of traps grows, as illustrated in the inset of Fig. 
3. A transition from insulating to conducting behavior is 
thus expected at a larger film thickness than the critical 
thickness for surface defect stabilization. With growing 
thickness, not only dopant levels tend to overlap as in 
a degenerate semiconductor, but the doping-level band 
is pushed towards the conduction band. For Z — 2, as 
for O vacancies, the physics of this transition is that of 
band overlap and disorder, since all dopant states are 
doubly occupied. If the mechanism involves Z = 1 de- 
fects, as in the hydroxylation case, the transition will be 
rather Mott-Anderson, as each dopant state is singly oc- 
cupied. The different phenomenologies could be used to 
ascertain on the mechanism. The surface potential dis- 
tribution from charged defects is consistent with a recent 
Kelvin probe force microscopy studjl^l 



IV. FINAL REMARKS AND CONCLUSIONS 

A recent AFM study of LAO/STO has proposed the 
mechanism for conductivity switching 4 5 as the writing 
of surface charge^. Applying a biased tip to the surface 
alters the field across the LAO film which either increases 
or decreases the stability of vacancies (and hence a v ) de- 
pending on the sign of the bias. An implication of this 
observation is that the kinetics for these redox processes 
is accessible at room T as used in these experiments, not 
only the much higher T used for growth. 




FIG. 3: (Color online) Trapping potential V created by a sur- 
face O vacancy as seen by interface electrons versus distance 
within the interface plane x (x=0 is directly below the va- 
cancy) for film thickness d — a (deepest), 2a, 3a, 4a (shallow- 
est), with harmonic estimates of corresponding donor ground 
states (taking m e fi = 3m e ?2 ). Inset: Sketch of range and 
density of trapped states. 



The model proposed can also be used for p interfaces, 
holes, acceptor levels and surface oxidation processes. 
This would be the case for LAO grown on SrO termi- 
nated STO. It is less symmetric than it seems, however, 
since, in addition to different energetics and chemical po- 
tentials, the large conduction band offset at the interface 
(W in Fig. 1) favors the situation for electrons towards 
n interfaces much more than the much smaller valence 
band offset for holes and p-interfaces. For thin film ferro- 
electrics with outwards (inwards) polarization on metallic 
substrates, the important alignment becomes the accep- 
tor (donor) level with the metal fermi level. This could be 
behind the stability of switchable ultra thin fer roelec tric 
films under open circuit conditions (see e.g. ref! 24 * 34 j). 

We conclude that in LAO/STO, the onset of electro- 
static modulation doping is precluded by the thermody- 
namic creation of surface defects and thus carrier mobili- 
ties produced by this method will be much lower than at 
a pristine interface. Intrinsic systems will likely require 
changed growth conditions or modified materials with a 
higher vacancy free energy. 
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Appendix: Model parameters 

Here we describe the determination of the parameters, 
a, e and C used in the model in Fig. 2 of the main 
paper. The parameters used in the model were inde- 
pendently determined from first principles in appropriate 
L AO-based systems (see below), and then the model was 
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checked against DFT calculations of vacancy formation 
energies in the full LAO/STO system as a function of 
thickness (Fig. 2a)). When comparing the model with 
experiment (Fig. 2b)), we account for inaccuracies of 
DFT and the ambiguity of the experimental chemical po- 
tential in the determination of each of these parameters. 

1. The vacancy- vacancy interaction term, a 

The defect-defect term was defined in Eq. 1 to include 
interactions other than electrostatic. Therefore to de- 
termine a we performed first principles calculations (see 
reP^ for the method) of the charge neutral defect, i.e. 
oxygen vacancies in 'bulk' LAO which include the dou- 
ble donor electrons. One oxygen vacancy was placed in 
a simulation cell of 1x1x8, 2x2x8 and 3x3x8 unit cells 
of LAO to approximate 2D arrays of vacancies of various 
area densities. From the difference in formation energy 
per vacancy between the three calculations, a was found 
to be 0.8 eV/(vac/f.u.) 2 , which was used in Fig 2a). a is 
formally defined as the interaction between vacancies at 
the film surface, however we believe this bulk value to be 
a good estimate. For the comparison with experiment, to 
account for any error associated with this determination 
we choose the range < a < 8 eV/(vac/f.u.) 2 in Fig. 2b) 
of the paper. 

2. The dielectric constant, e 

The dielectric constant, e, consists of lattice and elec- 
tronic contributions. For LAO, we take e=28 as the total 
for Fig 2a), as consistent with refP^. When compar- 
ing with experiment, we note the error and inconstancy 
of DFT calculations of e, and additionally the effect of 
strain as highlighted in ref.^. Due to this we choose the 
range 21 < e < 46 for Fig 2b). 

3. The surface/interface chemistry term, C 

C consists of three terms, 

C = E f ^ +» + 2(E cv -W) (A.l) 



From the electronic structure presented by Li et al., 
the last term is found to be approx. -1.2 eV. From first 
principles calculations of refp3 the formation energy of 
an isolated oxygen vacancy at the surface of LAO (in the 
absence of a field) with reference to oxygen in an iso- 
lated molecule (1/2S[0 2 ]), E*j^ =0 , is approx. 6.0 eV. 
We define the zero of chemical potential relative to this 
reference state, which is appropriate for the DFT com- 
parison and hence the value of C used in Fig. 2a) is taken 
as 4.8 eV. 

The DFT underestimation of the band gap requires 
corrections to both E® ^ =Q (see refP^l) and W for the 
comparison with experiment in Fig. 2b). At this point 
we note the difficulties and variation in first principles 
determination of formation energies of donor/acceptor 
states (see for example refP^). 

From refF^, the formation energy correction of a donor 
defect, required due to DFT band gap underestimation 
is simply: 

AE° f = ZAE C (A.2) 

where AE C is the change in conduction band edge be- 
tween LDA and experiment (or corrected DFT) . By com- 
paring the electronic structure presented in Li et al. and 
refPS, this correction could be as large as 1.0 eV. There- 
fore we take 6.0 eV < Ej^ < 7.0 eV. 

From experimental band alignment!^ and theoretical 
calculations determining the gap states^, the third term 
in Eq. (A.l) is approx. 2.0 eV (not 1.2 eV). Correcting 
for these DFT errors we take, 

4.0 eV + fx < C < 5.0 eV + fx (A.3) 

The chemical potential of oxygen in the growth con- 
ditions used in refP (T = 1073 K and p = 2.0 x 10~ 8 
atm) relative to the zero reference defined above is calcu- 
lated to be -1.9 eV assuming the environment acts as an 
ideal gas-like reservoir. The effect of post-annealing and 
cooling to room temperature and pressure is to shift the 
chemical potential towards zero. With these limits on 
the chemical potential and the inequality in Eq. (A.3), 
the range of C becomes 2.1 eV < C < 5.0 eV, as used in 
Fig. 2b). 
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